First order superconducting transition near a ferromagnetic quantum critical point 
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We address the issue of how triplet superconductivity emerges in an electronic system near a 
ferromagnetic quantum critical point (FQCP). Previous studies found that the superconducting 
transition is of second order, and T c is strongly reduced near the FQCP due to pair-breaking effects 
from thermal spin fluctuations. In contrast, we demonstrate that near the FQCP, the system avoids 
pair-breaking effects by undergoing a first order transition at a much larger T c . A second order 
superconducting transition emerges only at some distance from the FQCP. 

PACS numbers: PACS numbers: 74.25.-q 



Superconductivity near a magnetic instability is a topic 
of current interest in condensed-matter physics. Magnet- 
ically mediated pairing near an antiferromagnetic insta- 
bility is a candidate scenario for d— wave superconductiv- 
ity in the cuprates and heavy fermions compounds (for a 
recent review see (T)). The emergence of superconductiv- 
ity, mediated by the exchange of ferromagnetic spin fluc- 
tuations, is also expected near ferromagnetic transitions. 
Ferromagnetic exchange yields Cooper pairs with S = 1 
and therefore generally gives rise to triplet superconduc- 
tivity. This type of pairing was originally suggested by 
Anderson and Morel M for 3 He. In recent years, an in- 
tensive search has focused on superconductivity in com- 
pounds which can be tuned to a ferromagnetic quantum 
critical point (FQCP) by varying either pressure or chem- 
ical composition. Among the studied systems are MnSi, 
and the heavy fermion compound UGe2 (for an experi- 
mental review, see Ref. ||). 

The emergence of superconductivity in electronic sys- 
tems close to a ferromagnetic instability has recently 
been studied by three groups, who solved a linearized 
gap equation within the Eliashberg formalism Q |5[ ||. 
For both two- (2D) and three-dimensional (3D) systems, 
their analysis yielded a superconducting transition tem- 
perature, T l c (T stands for linearized), that substantially 
decreases as the system approaches criticality, eventually 
vanishing at the FQCP. 

The physical origin of the decrease in T\ near the 
FQCP lies in the presence of thermal spin fluctuations 
which behave like magnetic impurities whose scattering 
potential diverges as the critical point is approached ||, 
M . This behavior is reflected in the fermionic self-energy 
in the normal state, cx iT^ m sign (w n )xL(^m — 

ui n )i where Xl(u) = J d D ~ 1 qx{q,u) is the "local" spin 
susceptibility. Since Xl{^ = 0) diverges at the FQCP 
for D < 3 (assuming x(<?,0) = Xo/(£~ 2 + Q 2 ) with 
£ — > oo), the dominant contribution to E(w rl ) comes from 
the n = m term in the frequency sum, i.e., from classi- 
cal, thermal spin fluctuations. These fluctuations scat- 
ter with a finite momentum and zero frequency transfer; 
hence their similarity with magnetic impurities. This 
leads to = i-fsignui with 7 oc T£ 3 ~ D for system 



dimension D < 3 and 7 oc Tlog£ for D = 3. 

The analogy with magnetic impurities extends to the 
pairing problem in which thermal spin fluctuations close 
to the FQCP tend to break Cooper pairs and hence lower 
the temperature of the superconducting transition 0. 
This strong pair breaking effect is reflected in the gap 
equation where the divergent contributions to the self- 
energy E(w) and to the pairing vertex are not cancelled 
out in the equation for the pairing gap A(w) B (the ratio 
of divergent terms is 3 to 1 in our case). Simple estimates 
show that the linearized gap equation does not have a 
solution above 7 ~ T c q where T c q is the transition tem- 
perature in the absence of thermal fluctuations. Hence, 
when the FQCP is approached, T l c vanishes as T l c cx £ D ~ 3 
in D < 3, and T l c cx l/log£ for D = 3. Numerical so- 
lutions of the linearized Eliashberg equations near a fer- 
romagnetic instability demonstrate precisely this kind of 
behavior - T l c falls off when the FQCP is approached - 
more rapidly in 2D (§ than in 3D §, | (Refs. § § 
obtained a small finite T l c at criticality by using a self- 
consistent approach that goes beyond Eliashberg theory 
for the self-energy, but still neglects vertex corrections). 

In this letter, we argue that the actual behavior of the 
system is different from that discussed in Refs. ||, [|. 
Specifically, we show that close to a FQCP, supercon- 
ductivity emerges via a first order phase transition at 
T c ~ T c q. The much smaller T\ previously obtained by 
solving the linearized gap equation is just the end point 
of the temperature hysteresis loop, at which the normal 
state becomes unstable. The first indication that the 
pairing problem near the FQCP is unconventional comes 
from the observation that at T = 0, dangerous thermal 
fluctuations are absent, and hence the pairing gap A(a>) 
should generally be of the order of T c q. Explicit calcu- 
lations confirm this (see below). The second indication 
is that in the presence of a large gap, the spectrum of 
ferromagnetic spin fluctuations changes due to feedback 
effects from the pairing. For the ABM phase of triplet su- 
perconductivity, which we consider in the following, this 
feedback is different for Xzz and x± (assuming that 
the spin of the Cooper pair, S cp , lies in the xy-plane). 
At the FQCP, massless excitations survive in x±i but 



2 



not in Xzz which in the presence of a pairing gap de- 
scribes massive longitudinal spin fluctuations. Due to 
this distinction between x± arL d Xzz-, the ratio of the 
divergent terms in the self-energy and the pairing ver- 
tex is [2x4(0) + X \M\ I [2X^(0) - x\M\ ■ H X± and 
Xz Z behaved identically, this ratio would be 3, and pair- 
breaking effects of thermal fluctuations would be crucial. 
When only x± (0) diverges, the above ratio is 1 and the di- 
vergent terms from the fermionic self-energy and the pair- 
ing vertex cancel out in the gap equation. This in turn 
implies that the superconducting state with a large gap 
remains stable well above T°.. This behavior, however, 
cannot extend to small A(w) since then Xzz(0) cannot 
be neglected, the above ratio of divergent terms becomes 
3 and no cancellation occurs. We therefore expect that at 
the FQCP, the solution with a finite gap should survive 
up to the end point at Tf ~ T c0 > T l c ('nT stands for 
the solution of the nonlinear equation) where it becomes 
unstable. In other words, over some range of T, both the 
normal state and the state with a large gap are locally 
stable. This is a classic scenario for a first order transi- 
tion. Note, that this behavior is very different from that 
for singlet pairing, mediated by antiferromagnetic spin 
fluctuations, where the divergent contributions to H(u>) 
and to the pairing vertex cancel each other even in the 
linearized gap equation. As a result, the divergence of 
the thermal self-energy at criticality does not affect T l c 
H which saturates at a finite value for £ = oo. 

In the remainder of this paper we compute T™ 1 from 
the full set of nonlinear Eliashberg equations and show 
that it saturates at a finite value at the FQCP. We will 
not attempt to compute the actual T c (this would require 
the analysis of the condensation energy) . However, since 
T l c vanishes at the FQCP and T™ 1 stays finite, T c near 
criticality should be comparable to T" . 

Our starting point is the spin-fermion model, which de- 
scribes the interaction of low-energy fermions with their 
own spin degrees of freedom, S q , whose propagator is 
peaked at q = 0. The same model was used in earlier 
studies |, | 1 1, [l0[ . We assume that T c0 is much less 
than the Fermi energy Ep, implying that the pairing in- 
stability involves only fermions near the Fermi surface. 
The model is described by the Hamiltonian 

H = ^v F (k-k F )4 Q c k , Q + ^x" 1 (q)S q S_ q 

k,Q q 

+9 C k+q,a a °<-3 C KP ' S -q • M 

q,k,a,/3 

where the spin-fermion coupling, g, the Fermi velocity, 
vf (we assume a circular Fermi surface), and the static 
spin propagator x(q, 0) are input parameters. While 
the upper energy cutoff is in general also an input pa- 
rameter, our results are cutoff independent for the pair- 
ing problem considered here. The dynamical part of 
X(q, fi m ) = Xo/[q 2 + + n(g, Q m )] arises from the in- 
teraction with the low-energy fermions and is explicitly 
calculated. While we restrict our consideration to D = 2, 
our conclusions are also valid for 3-D systems. 



We assume that the static xil, 0) has the conventional 
lorentzian form with a weakly temperature dependent 
£. This form of x{l, 0) was recently questioned [jfl| [ij] 
since far away from criticality, the static spin susceptibil- 
ity possesses singular low-energy Fermi liquid corrections 
that give rise to a universal \q\ dependence of x{o.i 0), and 
a T dependence of £ . It is unclear, however, whether 
these singular corrections survive in the quantum criti- 
cal regime, so we restrict with the conventional form of 
x(q, 0) without further justifications. 

Near the critical point, a conventional perturbation 
theory in the spin-fermion coupling (for which £ = 
II = is the point of departure) holds in powers of 
A = .g 2 xo/(47rf_F^ _1 ), i.e., the quantum-critical region 
falls into the strong coupling limit. An approach for 
dealing with a strong coupling problem is the Eliashberg 
theory fl3|| . Its validity requires certain conditions to be 
met. We proceed assuming that the Eliashberg theory is 
valid, and then discuss what restrictions are necessary. 

We first consider the situation right at the FQCP 
where = 0. In the normal state, the dynamical 
part of the spin polarization operator, II(g, Q m ), is in- 
dependent of the fermionic self-energy, E(g,w n ) (but not 
vice versa), as the essential momenta for H(q,Q, m ) are 
those with vpq ^> £l m , S. As a result, IT(<7, O m ) has the 
same form as for free fermions, i.e., for |f2 m | <C vpq, 
U(q,n m ) = F(Q m )/(v F q) where F(fi m ) = ak%\Q m \, 
a = g 2 xo/(^Ep), and Ep — kpvp/2. At the same 
time, is determined by II and given by £(u> m ) = 

where cj = (3V3/A)a 2 E F . The non-Fermi 
liquid, w 2 / 3 -dependence of the self-energy in 2D is due to 
the divergence of the perturbation theory at the FQCP. 
This form was earlier obtained in Ref. |14|j . 

In the superconducting state, the equations for two 
components of F(£l) (F zz and F xx — F yy ) are coupled 
to the equation for the pairing gap A(w). Similar to the 
normal state, none of these quantities explicitly depends 
on S. As a result, one needs to self-consistently solve a 
set of three coupled equations for the two components 
of F and A. The derivation of the Eliashberg equations 
is quite straightforward and will not be presented here. 
With the Cooper pair spin lying in the xy-plane, the 
coupled equations for A(w), F ZZ (Q) = F_ and F xx = 
Fyy = F + at the FQCP have the form 

A K) = ^ T£ 



2u r , 



A(w m ) A (uj n ) 



[F+(T,A, 



1/3 



A(w m ) A (uj n ) 



[F_(T,A,< 



,1/3 



'm ) 



(2) 



and 



F±(T,A,w m -w„) = 

w n (w m ) ± A(w„)A(w m ) 



(3) 
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As anticipated, F + (T, A,0) = vanishes implying that 
the corresponding susceptibilities, Xxx and Xyy, describe 
massless modes. The vanishing of F + (T, A,0), how- 
ever, is not dangerous as it is compensated by the si- 
multaneous vanishing of the numerator in Eq.(^|). Ex- 
actly the same cancellation of divergences occurs in the 
gap equation for d-wave pairing due to antiferromag- 
netic spin fluctuations. The vanishing of F_ (T, A, 0) 
would be dangerous, but for finite A, F_(T, A,0) = 

2irTJ2 n A 2 (w„)/V w ™ + A2 ( w ») is also nnite > Le - the 
longitudinal spin excitation described by xZzi 1 ^ = 0) oc 
g 2 +F_ (T, A, 0)/(vFq) is massive. In contrast, for the lin- 
earized gap equation A is vanishingly small, F_ (T, A, 0) 
vanishes, and the r.h.s of Eq.(||) diverges. Due to this 
divergence, the linearized gap equation does not have a 
solution down to T — 0. Note, that the only energy scale 
in the Eliashberg equations is loq, which can be eliminated 
by rescaling both temperature and the gap in units of u> . 
The gap equation is then fully universal, which implies 
that the mass in \zz and the typical momentum q for the 
pairing problem are both of order ujq/vf- 

In Fig. |l|a we present the numerical solution for A(T) 
at the lowest Matsubara frequency, u> m = nT. As ex- 
pected for a first order phase transition, the gap changes 
discontinuously from a finite value to zero at T™ 1 ~ 
0.015cl>o- The inset shows that the discontinuous jump 
in A occurs for all Matsubara frequencies. 

We next study the situation at finite £ and verify 
whether the first order superconducting transition be- 
comes second order at some distance from the FQCP. 
Away from criticality, the equations for F±(T, A,u> m ) 
retain their form, but in the gap equation, the fac- 
tors F±(T,A,w m ) are replaced by F±(T, A, uj m )/I 3 (fl±), 
where I{(3±) by itself depends on F±(T, A, u m ) through 



3\/3 
~2tT 



dx- 



1 + (3±x + , 



(4) 



Here (3± = b/[F±(T, A, il m )/uj Q ] 2 / 3 , and b = 
(8/3y/3) 2 ^ 3 (akp^)~ 2 measures the deviation from the 
FQCP (b = at the FQCP). The integral can be ex- 
pressed in elementary functions. One can easily verify 
that for finite £, the gap equation does not contain any 
divergence, even at infinitesimally small A, and hence 
both T™ ( and T c are nonzero. 

In Fig. |l|b we plot A(T, inT) for 6 = 2. We clearly 
see that the temperature dependence of the gap is now 
continuous, in marked distinction to Fig. [pa. The inset 
shows that the continuous evolution of A holds for all 
Matsubara frequencies. This result implies that for large 
enough b > b c , the transition is of second order. To lo- 
cate the tricritical point b = b c , we plot in Fig. |l|c the 
magnitude of the jump of A(T, inT) at T™ 1 as a func- 
tion of b. We see that the gap discontinuity disappears 
at b c w 1.2. For completeness, in Fig. |l|d we plot the 
zero temperature gap versus b. We see that it changes 
gradually, without a singularity at b c . 

The phase diagram that emerges from our studies is 
presented in Fig. 0. The actual superconducting tran- 
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FIG. 1: (a) Temperature dependence of A(T, u> m ) at the low- 
est Matsubara frequency oj m — nT for b = 0. The lines are 
a guide for the eye. A, T and ui m are in units of coo (see 
text). The discontinuity of A(T) at 0.015 indicates a first 
order transition. The inset shows A(a>) versus frequency at 
several T. (b) Same away from the FQCP, for b — 2. Now 
the transition is continuous, (c) The magnitude of the jump 
of A(T, inT) at the instability temperature versus b. The line 
is a guide to the eye. (d) A(T, inT)) at the lowest T versus b. 



sition temperature lies between T l c and T™ z and hence 
remains finite at b = 0. For 0.9 < b < b c , the jump in A 
at T™ 1 , and consequently the difference between T l c and 
T™ ( is small but finite (see Fig. |l|c). The insert shows 
the reduction of T c at large b due to the decrease in the 
effective coupling. 

The functional form of T l c (b) obtained above agrees 
qualitatively with that found earlier by Monthoux and 
Lonzarich Q]. Their result for the maximum T l c is 
T lmax ~ 0-016 lj b 3 ' 2 . Substituting b - 0.8, where T l c is 
maximal in our case, we obtain T l c max ~ O.Ollwo, con- 
sistent with our result T l c max ~ 0.013ljo- 

Finally, the validity of the Eliashberg treatment re- 
quires that three conditions be satisfied. First, typical 
bosonic momenta qs should be much larger than typical 
fermionic \k— fopl, i.e., bosons should bc slow modes com- 
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FIG. 2: The phase diagram near the FQCP. In the near vicin- 
ity of the FQCP, the transition is of first order, away from the 
FQCP, to the right of b c , it is of second order. For the first 
order transition, T" 1 and T c are the instability temperatures 
for the solutions with a large and infinitesimally small gap, 
respectively. The actual first order transition temperature, 
T c , lies between T l c and T™ 1 . Inset: the reduction of T c at 
large b. 

pared to fermions, leading to £(£;, w) « A straight- 

forward analysis shows that the typical qs ~ akp, while 
the typical \k — kp\ ~ uJo/vp ~ a 2 kp. The Eliashbcrg 
theory is therefore valid if a <C 1, i.e., ojq <C Ep. This 
in turn implies that the physical behavior of the system 
is universally determined by only low energy excitations. 
Second, vertex corrections should be small. Generally, 
this is not possible for typical qb kp, as vertex correc- 
tions scale with £. However, for a -C 1, we only require 
the vertex for O -C vpq since the typical «f9b well ex- 
ceeds the typical Q ~ ujq. In this limit, vertex corrections 
are much smaller and only scale as log£. They are still 



non-negligible at the FQCP, where they change the pole 
in the spin susceptibility into a branch cut [jlo|| . However, 
we verified that, as in the antiferromagnetic case ||, this 
only leads to a small renormalization of the prefactors in 
the gap equation. Third, one should be able to neglect 
the momentum dependence of the pairing gap, while pre- 
serving the gap symmetry, A(nkp) = —A(—nkp). This 
approximation is again justified by a <C 1, as in this limit, 
the typical momentum transfers along the Fermi surface 
8k r*> qs ~ ockp are much smaller than kp. In this situ- 
ation, the momentum variation of the gap at typical 5k 
only introduces 0(1) corrections to the Eliashberg theory 
which can be safely neglected. Note in passing that 
the smallness of 8k <C kp makes our theory also appli- 
cable to real materials (in which a crystalline structure 
imposes additional constraints on the order parameter 
symmetry ) , as it allows one to consider the pairing 
problem in a local-momentum approximation, ignoring 
the peculiarities of the gap's momentum dependence. 

In summary, we showed that near a FQCP, spin fluc- 
tuation exchange gives rise to a strong first order transi- 
tion into a p— wave superconducting state. By choosing 
a first order transition, the system avoids divergent pair- 
breaking effects from thermal spin fluctuations. As a 
result, T c saturates at a nonzero value at criticality. The 
first order transition persists up to a finite distance from 
the FQCP, where it becomes second order. 
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